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We introduce a multi-species network model which describes the resistance fluctuations of a 
resistor in a non-equilibrium stationary state. More precisely, a thin resistor characterized by a l/f a 
resistance noise is described as a two-dimensional network made by different species of elementary 
resistors. The resistor species are distinguished by their resistances and by their energies associated 
with thermally activated processes of breaking and recovery. Depending on the external conditions, 
stationary states of the network can arise as a result of the competition between these processes. The 
properties of the network are studied as a function of the temperature by Monte Carlo simulations 
carried out in the temperature range 300 -j- 800 K. At low temperatures, the resistance fluctuations 
display long-term correlations expressed by a power-law behavior of the auto-correlation function 
and by a value ~ 1 of the a-exponent of the spectral density. On the contrary, at high temperatures 
the resistance fluctuations exhibit a finite and progressively smaller correlation time associated with 
a non-exponential decay of correlations and with a value of the a-exponent smaller than one. This 
temperature dependence of the a coefficient reproduces qualitatively well the experimental findings. 

PACS numbers: 05.40.-a, 72.70.+m, 61.43.-j, 64.60.Ak 
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I. INTRODUCTION 

Since a long time the analysis of the resistance fluctua- 
tions has turned out to be a very powerful tool for probing 
various condensed matter systems uiM an d, in partic- 
ular, disordered materials [3, 0, S H S] like conductor- 
insulator composites [|| , granular systems Q , porous [1] 
or amorphous materials [ToL fllj and organic conduct- 
ing blends Therefore, many experimental and the- 
oretical investigations have been devoted to study the 
properties of the resistance noise as a function of tem- 
perature, bias strength and of the main material param- 
eters [3, H !, 0, S 1 M EL Ei 0, 0, EE El Ep- 

Only few studies have considered the behavior of dis- 
ordered materials over the full range of bias values, by 
analyzing under a unified approach the linear regime 
and the nonlinear regime up to the threshold for elec- 
trical breakdown [f| Ell El Hll, even though this kind 
of studies can provide significant insights about basic 
properties of the response of nonequilibrium systems 

EEiElEIIMIIlIilil. 

In previous works we have introduced [24j and de- 
veloped the Stationary and Biased Resistor Network 
(SBRN) model [19, 20, 25, 26] which is based on a resistor 
network approach [3|, [j, Q and which allows the study of 
the electrical conduction and noise of disordered materi- 
als over the full range of bias values. This model describes 
a disordered conducting material as a resistor network 
whose structure is determined by a stochastic competi- 
tion between a biased percolation [2?], HI] of "broken" 
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resistors (high resistivity resistors) and a biased recovery 
process (which restores the "regular" resistors with "nor- 
mal" resistivity) [HHi]. The SBRN model provides a 
good description of many features associated with the re- 
sistance fluctuations in none quili brium stationary states 
of disordered resistors Irj Tl9l I20I l25l |29| , with resis- 
tor trimming processes [301 ] and with agglomeration and 
thermal instability phenomena in metallic films [3ll . 1321 ] . 
Further successful applications of the model concern the 
study of electrical breakdown phenomena in composite 
and granular materials [6J, [l9|, l2fj, |25| . including electro- 
migration damage of metallic lines [l6t [26l [33] ] . 

However, this model only applies to systems character- 
ized by a Lorentzian noise i.e. by a power spectral 
density of the resistance fluctuations scaling as 1/f 2 at 
high frequencies and becoming flat at frequencies below 
the so called corner value f c . The existence of a cor- 
ner frequency in the spectrum corresponds in the time 
domain to the existence of a well defined characteristic 
time, r (correlation time), associated with the decay of 
the correlations in the resistance fluctuations In the 
particular case of Lorentzian noise this decay is exponen- 
tial, i.e. the auto-correlation function of the resistance 
fluctuations is given by [l[: 

C SR (t) = (SR(t)5R(t + r)) = ((AR) 2 } exp(-t/r) (1) 

where ((AR) 2 ) is the variance of the resistance fluctu- 
ations. On the other hand, it is well known [l[ that 
at low frequencies many condensed matter systems dis- 
play 1/f resistance noise, i.e. a spectral density scal- 
ing over several frequency decades as 1/ f a , with a ~ 1. 
The amazing presence of 1/f noise in a large variety of 
natural phenomena represents a puzzling problem not 
yet fully solved. Therefore, a huge quantity of mod- 
els providing signals with 1/f noise have been proposed 
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particular, the ubiquity of the 1/f noise has given rise 
to many attempts to explain it in terms of a univer- 
sal law. One simple way to obtain a 1// spectrum is 
by superimposing a large number of Lorentzian spectra 
with an appropriate distribution of the correlation times: 
D(t) oc 1/r, distributed over many orders of magnitude 
of t values In some cases, the distribution of the 
correlation times of the elementary processes contribut- 
ing to the 1/f spectrum can be derived from the dis- 
tribution of some variable on which the correlation times 
themselves depend, as in the case of the pioneering works 
of Mc Whorter [l| and Dutta and Horn [l[ . In particu- 
lar, these authors proposed [l[ that the origin of the 1/f 
noise can be attributed to a thermally activated expres- 
sion of the correlation times, associated with a broad 
distribution of the corresponding activation energies, an 
assumption physically plausible for many systems. On 
the other hand, many other important contributions to 
the understanding of the 1/f noise have been advanced 
in the last twenty years, which have shown that the pres- 
ence of a 1/f spectrum can also arise from other basic 
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reasons 

The dynamical random network model 
spin-glasses models 0, [H{ provide a good explanation 
of the 1/f noise in conducting random magnetic mate- 
rials. Dissipative self-organised criticality (SOC) mod- 
els, started from the famous work of Bak, Tang and 
Wiesenfeld [35| , clarify the origin of 1// spectra in certain 
dissipative dynamical systems naturally evolving into a 
critical state. Avalanche models 13611. clustering models 
[37l | and percolative models @, UTIH, HE E3| represent 
other relevant classes of theoretical approaches explain- 
ing the appearence of 1/f noise in a variety of systems. 
Thus, the conclusion, now largely accepted in the liter- 
ature [H, 0, E(|, i s tnat a unique, universal origin of 1/f 
does not exist, even though classes of systems can share 
a common basic origin of 1// noise. 

The aim of this paper is to present a new model which 
includes and extends to systems characterized by l/f a 
noise the main feature of the SBRN model, i.e. the 
stochastic competition between thermally activated and 
biased processes of breakingand recovery of the elemen- 
tary resistors of a network pj| . Actually, this feature has 
turned out to be essential for the success of the SBRN 
model in the description of the wide phcnomclogy associ- 
ated with the resistance fluctuations of disordered resis- 
tors, in both nonequilibrium stationary states and non- 
stationary failure states 0, EE, US [H, H EH E2] . To this 
purpose, here we study a network made by several species 
of resistors, where each species is characterized by a resis- 
tance value and by a pair of activation energies associated 
with the breaking and the recovery processes. For this 
reason we call this model "multi-species network" (MSN) 
model. The states, either stationary or nonstationary, of 
the MSN result from the stochastic competition between 
these processes of breaking and recovery involving the el- 
ementary resistors. In analogy with the Dutta and Horn 



model [H, we take the breaking and recovery activation 
energies distributed in a broad range of values. A choice 
physically justified in disordered systems, like granular 
and amorphous materials, composites, etc. where the 
orientational disorder present inside these materials can 
give rise to different energy barriers for the electron flow 
along the different conducting paths [l], 0, 0, [13] • As a 
result, the resistance fluctuations of the MSN exhibit a 
1/ f a spectrum, where the value of a depends on the ex- 
ternal conditions (temperature and electrical bias). In 
particular, in this paper, we discuss the states of a MSN 
as a function of the temperature and in the linear regime 
of the external bias, when Joule heating effects are neg- 
ligible. 

The paper is organised as follows: in Sect. II we il- 
lustrate the MSN model, while in Sect. Ill we report 
the results by focusing in particular on the role played 
by the temperature on the fluctuations properties of the 
multi-species network. Finally, in Sect. IV we draw the 
conclusions of this study. 



II. MODEL 

According to the SBRN model [IE HE] a conducting 
thin film with granular structure is described as a two- 
dimensional square-lattice resistor network [E 0|. For 
simplicity here we study a network N x N (where N 
is the linear size of the network), even if networks with 
a different geometry can be considered 26]. The net- 
work is biased by an external constant current, /, applied 
through perfectly conducting bars placed at the left and 
right hand sides and it lies on an insulating substrate at a 
given temperature T, which acts as a thermal bath. Each 
resistor can be in two different states: i) regular, corre- 
sponding to a resistance r n = ro[l + otT(T n — T)], where 
(Xt is the temperature coefficient of the resistance and 
T n is the local temperature and ii) broken, correspond- 
ing to an effectively "infinite" resistance, tqp — 10 9 r„. 
Resistors in the broken state will be called defects. The 
temperature T n is expressed as [28| : 



T n = T + A 
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where D = 3/4, the sum is performed over the N neig 
nearest neighbors of the n-i/i resistor and i n is the cur- 
rent flowing in it. The above expression of T n takes into 
account the Joule heating of the n-th resistor [27I. |28| . |44T | 
and the thermal exchanges with its neighbors [28]. The 
importance of both these effects is controlled by the pa- 
rameter A, thermal resistance, which describes the heat 
coupling of each resistor with the substrate [IE, E3, EH • 
We note that by adopting Eq. $2$ we are assuming for sim- 
plicity an instantaneous thermalization of each resistor 
and we are neglecting time dependent effects discussed 
in Ref. 0. 

In the initial state of the network all the resistors are 
taken to be identical: r n = tq. The SBRM model then 
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assumes that the evolution of the network is determined 
by the competition between two biased stochastic pro- 
cesses: one of breaking and the other of recovery. The 
former process consists in the transition of an elementary 
resistor from the regular to the broken state and it occurs 
with a thermally activated probability [2^ |: 



of resistance values: ro,, £ [r. 



mm : ' max 



], while other op- 



W Dn = cxp{-E D /k B T n ) 



(3) 



the latter process consists in the reverse transition and it 
occurs with probability [24[ : 



W Rn = cxp{-E R /k B T n ) 



(4) 



where Ed and Er are the activation energies of the two 
processes and k B the Boltzmann constant. The time evo- 
lution of the network is then obtained by Monte Carlo 
simulations which update the network resistance after a 
sweep of breaking and recovery processes, according to 
an iterative procedure detailed in Ref. [l9j. The se- 
quence of successive network configurations provides a 
resistance signal, R(t), after an appropriate calibration of 
the time scale. Then, depending on the stress conditions 
(J and T) and on the network parameters (size, activa- 
tion energies, ro, ar and A), the network either reaches a 
steady state or undergoes an irreversible electrical failure 
[l9l I20I l25l . H(| . This latter possibility is associated with 
the achievement of the percolation threshold, p c , for the 
fraction of broken resistors [Sj]. Therefore, for a given net- 
work at a given temperature, a threshold current value, 
Ib , exists above which electrical breakdown occurs . 
For values of the applied current below this threshold, 
the steady state of the network is characterized by fluc- 
tuations of the fraction of broken resistors, Sp, and of the 
resistance, SR, around their respective average values (p) 
and (R). 

The network described above, as provided by the 
SBRN model, is made by a single species of resistors, 
all characterized by the same values of ro, ar, Ed and 
Er. Thus, we can speak of a "single-species network". 
The disorder inside this network only arises from the 
fact that some resistors can be in the broken state. As 
a consequence, current and temperature are not homo- 
geneously distributed [2(| H?| , thus implying inhomoge- 
neous probabilities of breaking and recovery and, through 
the temperature coefficient ax, different resistance val- 
ues of the regular resistors. Furthermore, by analogy 
with generation-recombination noise in a two-levels sys- 
tem, this model describes a system characterized by a 
single time scale and thus by Lorentzian noise. To over- 
come the above limitations here we consider a network 
made by N spec species of resistors, i.e. a "multi-species 
network" (MSN). Each species is characterized by differ- 
ent values of the parameters: ro,i, 0-T,h £x>,i and Er^ 
with i = 1, ...N spec . Then, an "identity" (an index i) is 
attributed at random to each resistor. In the initial state 
of the network the species are taken to be present with 
the same number of resistors. As a reasonable choice, we 
have taken ro,i uniformly distributed inside a given range 



tions are as well reasonable. We limit the present study to 
the linear regime in the external bias, when Joule heating 
effects are negligible, i.e. T n = T, Vn and consequently 
r n j = ro.i, Vn and Vi. Alternatively, we can say that we 
study a system characterized by negligibly small values 
of the thermal resistance and of the resistance temper- 
ature coefficients: A — and ar.% = 0, Vi. We stress 
that the linear regime condition corresponds to uniform 
probabilities of breaking and recovery for the resistors 
belonging to each given species [H, SBJ- Therefore, the 
results that we will rep ort in Sect. Ill directly generalize 
the results of Ref. [45| to systems characterized by 1// Q 
resistance noise. The limitation to the linear regime of 
currents has been adopted here for simplicity and also be- 
cause we aim to focus on the temperature dependence of 
the correlation properties of the resistance fluctuations. 
However, we underline that there is no special difficulty 
in extending this study to the nonlinear regime in the 
applied bias. In this respect, it must be noted that we 
expect that all the other significant features associated 
with the SBRN model (existence of a threshold current 
for breakdown [l9|, dependence of (R) and {(AR) 2 ) 
on the bias [lj|, times to failure distributio n |2ql . non- 
Gaussianity of the resistance fluctuations [2fll25l|. role of 
the network geometry (26j . etc.) remain unchanged for a 
MSN network. 

The values of the activation energies Ed,i and Er^ 
for the different resistor species have been chosen as fol- 
lows. In Ref. [45| it was proved that for a single-species 
network, when the average fraction of defects (p) is suf- 
ficiently far from the percolation threshold (p c = 0.5 for 
aJVxJV square-lattice network [f|), the correlation time 
of the resistance fluctuations is approximated by: 
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Wr = Wd 

(1 - Wr) (p) 



(5) 



where Wd and Wr are the thermal activated proba- 
bilities defined by the Eqs. ©-(HI) and calculated for 
T n = T. In other terms, Eq. ^ implies that at a given 
temperature, both (p) and r only depend on Ed and Er. 
This equation can thus be inverted to determine Ed and 
Er as functions of (p) and r. These considerations can 
be generalized to the case of a multi-species network by 
giving the criterion for a convenient choice of the activa- 
tion energies. To this purpose, first we choose a reference 
value of the substrate temperature, T re f, and a set of 
values of the average defect fraction, (p,;), corresponding 
to the different resistor species (where pi = N{,rok,i/N t ot,i 
with i = 1, N spec ). Second, by considering the correla- 
tion times, Ti, associated with the fluctuations [461 ] of Pi, 
we choose the values of the activation energies Ed,% and 
Er^ in such a way to obtain a logarithmic distribution 
of Ti inside a sufficiently wide interval: £ [T m i n ,T max \. 
The different Ti are then attributed to the different re- 
sistor species by adopting the criterion that increasing 
values of Tj are paired with increasing values of ro,i- Al- 
ternative choices are of course possible. 
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All the results reported here are obtained by applying a 
bias current I = 1.0 (A) to a network of sizes 75 x 75 made 
by N spec = 15 resistor species and by choosing r min = 0.5 
£1 and r max = 1.5 Q. Furthermore we take: T re f = 300 
K, (p^ ks 0.25 Vi, T min w 2 and T max w 5 x 10 5 (where 
times are expressed in units of iterative steps). We un- 
derline that the choice made above for the values of T re f, 
(pi), Train and T max , determines once for all (i.e. for all 
the temperatures) the set of activation energies Ed,* and 
Ejh. In the present case Ed i and Er^ are in the range 
58 -j- 375 meV and 37 -=- 346 meV, respectively. In fact, 
the values of Ed,i and Er^, together with ro,, and N spec , 
determine the intrinsic properties of the material under 
test. Therefore, when performing a simulation at temper- 
atures different from T re f, the values of r, recalculated 
from Eq. ([5]) are the input parameters of the simulation, 
while the values of (pi) only serve as a check of consis- 
tency with those obtained directly from the output of the 
simulation. As a general trend, by increasing the tem- 
perature the range of values of t; becomes more narrow 
and the values of (p^ increase progressively. 

The time evolution of the network is then obtained 
by Monte Carlo simulations according to the iterative 
procedure described in details in Ref. The auto- 

correlation functions and the power spectral densities of 
the resistance fluctuations are calculated by analyzing 
stationary signals R(t) consisting of 2 x 10 6 records. 



III. RESULTS 

The resistance evolution of a multi-species network cal- 
culated at a temperature of 300 K is reported in Fig. 1. 
The inset displays a small part of the same evolution over 
an enlarged time scale. The co-existence of different char- 
acteristic time scales in the R(t) signal is clearly shown 
in Fig. 1, where the long relaxation time associated with 
the achievement of the steady state, r re i , co-exists with 
the shorter times characterizing the resistance fluctua- 
tions displayed in the inset. For comparison, we report 
in Fig. 2 the time evolution of the resistance of a single- 
species network obtained at T = 300 K by the SBRN 
model. In this case, the values of the activation energies 
(Ed = 350 meV and Er — 310 meV) are chosen to give 
a relaxation time comparable with that of the signal in 
Fig. 1. The behavior of the resistance is now determined 
by a single time scale (r « T re {) and thus it is less noisy 
and essentially flat over time scales shorter than r. This 
is emphasized by the inset of Fig. 2 where the stochas- 
tic signal resembles that of a few level system (it must 
be noted that the vertical scale of the inset in Fig. 2 is 
significantly enhanced with respect to that of Fig. 1). 

Figure 3 reports the auto-correlation functions of the 
resistance fluctuations on a log- log plot. More precisely, 
the auto-correlation function corresponding to the multi- 
species network (MSN model) is shown by black trian- 
gles while that corresponding to the single-species net- 
work (SBRN model) is reported by black squares. The 
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FIG. 1: Resistance evolution of a multi-species network (MSN 
model) calculated at 300 K. The resistance is espressed in 
Ohm and the time in iterative steps. The inset highlights the 
resistance fluctuations on an enlarged time scale. In particu- 
lar, the time units are divided by a factor 10 . 
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FIG. 2: Resistance evolution of a single-species network 
(SBRN model) calculated at 300 K. The resistance is 
espressed in Ohm and the time in iterative steps. The in- 
set displays the resistance fluctuations on the same enlarged 
time scale of the inset in Fig. 1. (the vertical scales of the 
insets in the two figures are different). 

solid and short dashed lines represent the best-fits to 
the two Csr functions carried out respectively with a 
power-law and with an exponential law. The fitting pro- 
cedure confirms the exponential decay of the correlations 
in the resistance fluctuations of the single-species net- 
work, described by Eq. ([T]), and points out the existence 
of long-term correlations in the resistance fluctuations of 
the multi-species-network, characterized by a power-law 
decay of the auto-correlation function: 

C iB (t)~t-T (6) 

with < 7 < 1. In particular, for the case of Fig. 
3 we have found a value 7 = 0.22 for the correlation 
exponent. We stress that the above expression of the 
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FIG. 3: Auto-correlation functions of the resistance fluctu- 
ations calculated for a multi-species network (MSN model, 
black triangles) and for a single-species network (SBRN 
model, black squares). Both functions are obtained at 300 
K. The solid and short dashed grey lines show the best-fit 
respectively with a power-law of exponent 7 = 0.22 and with 
an exponential function with correlation time r = 5.8 x 10 4 . 
The time is expressed in iterative steps. 




Frequency (arb. units) 

FIG. 4: Power spectral density of the resistance fluctua- 
tions at 300 K calculated for the multi-species network (MSN 
model, solid line) and for the single-species network (SBRN 
model, dotted line). The grey solid line shows the best-fit to 
the MSN spectrum with a power-law of slope -0.94. The grey 
dashed curve represents the best-fit with a Lorentzian to the 
SBRN spectrum. 



auto-correlation function implies a divergence of the cor- 
relation time, as can be easily seen by considering the 
following general expression of the correlation time [47| : 



C SR (t) 

c SR (oy 



dt 



(7) 



SBRN model are displayed in Fig. 4. The grey solid 
line represents the best-fit with a power-law to the MSN 
spectrum. The grey dashed curve is the best-fit with a 
Lorentzian to the SBRN spectrum. The corner frequency 
of the Lorentzian, f c = 4.0 x 10~ 6 (arbitrary units), is 
consistent with the correlation time reported in Fig. 3 
and obtained by the best-fit of the corresponding auto- 
correlation function. Thus, Fig. 4 shows that at 300 K 
the resistance fluctuations of the multi-species network 
exhibit a power spectral density scaling as 1// Q over sev- 
eral decades of frequency, with a value a = 0.94. This 
result is a consequence of the envelope of the different 
time scales associated with the different resistor species. 

Now, we will analyze the effect of the temperature on 
the resistance fluctuations of a multi-species network. To 
this purpose, Fig. 5 displays the resistance evolutions 
calculated at increasing temperatures: T = 400 K (lower 
curve) and T = 600 K (upper curve). Already the qual- 
itative comparison of the R(t) signals in Figs. 1 and 
5, shows that at increasing temperature there is a sig- 
nificant growth of both the average resistance and the 
variance of the resistance fluctuations (we will discuss in 
detail the behavior of these quantities in the following) . 
Moreover, these figures point out a drastic reduction of 
the relaxation time at increasing temperature (more than 
one order of magnitude when T rises from 300 K to 400 
K). 
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The power spectral densities of the resistance fluctua- 
tions calculated at 300 K by the MSN model and by the 



FIG. 5: Resistance evolution of a multi-species network at 
400 K and 600 K. The resistance is espressed in Ohm and the 
time in iterative steps. 

The temperature is also found to affect the distribution 
of the resistance fluctuations by increasing its skewness, 
as shown in Fig. 6, which reports on a lin-log plot the 
probability density function (PDF, here denoted by (f>), 
of the distribution of the resistance fluctuations. A nor- 
malized representation has been adopted for convenience 
(here a is the root mean square deviation from the aver- 
age resistance). In Fig. 6, the different symbols represent 
the PDFs calculated in the steady state of the multi- 
species network at 300 K (full diamonds), 400 K (open 
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FIG. 6: Normalized probability density functions of the re- 
sistance fluctuations of a multi-species network calculated at 
300 K (full diamonds), 400 K (open circles), 500 K (open 
up-triangles) and 600 K (full down-triangles); a is the root- 
mean-square deviation from the average resistance value. The 
solid black curve shows the normalized Gaussian distribution. 

circles), 500 K (open up-triangles) and 600 K (full down- 
triangles). For comparison, the normalized Gaussian dis- 
tribution has been also reported as a solid line. The figure 
shows a significant non-Gaussianity of the resistance fluc- 
tuations, which becomes stronger at high temperatures, 
meaning that the system is approaching failure condi- 
tions. The distribution of the resistance fluctuations and 
the role played on the non-Gaussianity by the size, shape 
and disorder of the network has been investigated in Rcfs. 
(20l |25| , where the link with the universal distribution of 
fluctuations of Bramwell, Holdsworth and Pinton [48|, H^] 
has been discussed. 

Figure 7 reports the auto-correlation functions of the 
resistance fluctuations of a multi-species network calcu- 
lated at 400 K and 600 K. For comparison, the CgR func- 
tion obtained at 300 K (already shown in Fig. 3) has 
been drawn again in Fig. 7 with the solid grey line still 
giving the best-fit with a power-law. The dashed grey 
curves instead represent the best-fit to the Csr functions 
at 400 K and 600 K with the expression: 

C SR (t) = C t-' l eM-t/u) (8) 

The values of the best-fit parameters are: Co = 1.13, 
h = 0.30 and u = 1.42 x 10 4 for T — 400 K and 
C = 0.965, h' = 0.46 and u' = 5.57 x 10 2 for T = 600 
K. As pointed out by Fig. 7, the fit to C S r with Eq. © 
is very satisfying. We thus conclude that at high tem- 
peratures the auto-correlation function of the resistance 
fluctuations of the multi-species network is well described 
by a power-law with an exponential cut-off. Such kinds 
of mixed decays of the correlations, non-exponential and 
non power-law, are often found in the transition of a 
complex system from short-term correlated to long-term 
correlated regimes (5(1 IHl[ . It should be noted that for 



u — > 00 Eq. ||5J) becomes a power-law while for h — > it 
describes an exponential decay. 

Actually, Fig. 7 makes evident a strong reduction of 
the correlation time of the fluctuations when the temper- 
ature increases. This result can be easily understood in 
terms of Eq. ([5]) , by considering the thermally activated 
expressions of the breaking and recovery probabilities: in 
fact the increase of temperature implies a significant re- 
duction of the ratio T m i n /T max . Therefore the interval 
[Tmini T m ax] where the Tj are distributed becomes more 
and more narrow, destroying the power-law decay of cor- 
relations and giving rise to an exponential tail. It must 
be noted that, for similar reasons, when the temperature 
decreases below the reference value, T < T re f, the ratio 
Tmin/Tmax strongly increases and the correlations in the 
resistance fluctuations keep their power-law decay over 
wider time scales (for example, when T — 200 K the ra- 
tio T m i n /T max becomes 7.5 x 10 7 ). T re f is thus closely 
related to the transition temperature, T* , from a long- 
term correlated behavior i.e. a behavior characterized by 
a divergent correlation time and occurring for T < T* , 
to a behavior characterized by a finite correlation time 
and occurring for T > T* . 

The correlation time of the resistance fluctuations can 
be directly estimated at T > T re f by considering Eq. 
(7) and making use of Eq. (jHJ). We obtain the following 
expression of r in terms of the best-fit parameters of the 
auto-correlation function: 

t = u l - h T{l~h) (9) 

where Y is the Gamma function. The values of r calcu- 
lated in such a way are reported in Fig. 8 as a function 
of the temperature. The figure shows a strong increase 
of r at decreasing temperatures, when T approches T re f, 
in fact for this value of temperature r is expected to 
diverge consistently with the power-law behavior of the 
auto-correlation function. To highlight the dependence of 
r on the temperature, we report in Fig. 9 a log-log plot of 
the correlation time as a function of the difference T— T* , 
where the value of the transition temperature T* = 306 
K is determined by a best-fit to the data in Fig. 8 with 
the power-law: r ~ (T — T*)~ e . The fitting curve is dis- 
played in Fig. 9 by the dashed straight line. We find a 
value 9 = 2.7 for the exponent. Therefore Fig. 9 shows 
that the dependence of the correlation time on temper- 
ature is well described by a power of T — T* , where the 
transition temperature T* ks T re f. 

The power spectral densities of the resistance fluctu- 
ations calculated at the temperatures of 400 K and 600 
K are reported in Fig. 10. Here the grey lines represent 
the best-fit with power-laws of slopes -0.87 and -0.78, re- 
spectively. Thus we find that the power spectrum keeps 
the l// a form also at high temperatures. However, the 
value of the exponent a decreases to values well below 
unity when T > T* . This dependence of a on tempera- 
ture has been studied since a long time Pi [HI- A quali- 
tative description of this dependence is provided by the 
well known Dutta-Horn relation [l| . On the other hand, 
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FIG. 7: Auto-correlation functions of the resistance fluctua- 
tions of a multi-species network calculated at different tem- 
peratures. The solid grey curve shows the best-fit with a 
power-law to the auto-correlation function at 300 K (the same 
of Fig. 3). The dashed grey curves display the best-fit to the 
auto-correlation functions at 400 and 600 K with the func- 
tion: C(t) = Cot~ h exp[— t/u\ (see the text for the values of 
the fit parameters). The time is expressed in iterative steps. 
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FIG. 8: Correlation time of the resistance fluctuations as a 
function of the temperature. The time is expressed in iterative 
steps and the temperature in K. The dotted line connecting 
the diamonds is a guide to the eyes. 



FIG. 9: Correlation time of the resistance fluctuations as a 
function of the difference T — T* . The time is expressed in 
iterative steps and the temperature in K. The value of T* is 
reported in the figure. The dashed line shows the fit with a 
power-law of exponent 8 — 2.7. 
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FIG. 10: Power spectral density of the resistance fluctuations 
of the multi-species network at T = 400 K and T = 600 K. 
The grey lines show the best-fit with power-laws of slopes 
-0.87 and -0.78, respectively. 



many experiments have pointed out a strong dependence 
of the detailed behavior of a(T) on the particular mate- 
rial [l|, d, [aH- In this respect, it should be noted that a 
decrease of the a exponent from a ~ 1 and a « 0.8—0.5 is 
frequently observed in the experiments [l], 0, H3] at inter- 
mediate temperatures above 100-200 K, similarly to our 
case where this decrease occurs above T re f ss T* « 300 
K. It should be noted, however, that in our model T re f 
merely plays the role of an input parameter which can be 
adapted to fit the experiments. Furthermore, we stress 
that the monotonic decrease of a for T > T* is obtained 
here in the linear regime of currents, i.e. neglecting Joule 
heating effects. Actually, these effects, whose importance 



also depends on temperature can give rise to a more com- 
plicated behavior of a versus temperature. 

Now we consider the dependence on temperature of 
the average resistance and of the variance of the resis- 
tance fluctuations. We call Rq the average value of the 
resistance in the limit T — 0. This value represents the 
resistance of a network free of broken resistors but in any 
case disordered (due to the presence of different resis- 
tor species). Thus i?o not only depends on the network 
size and on the set of values {ro,i} with i = 1, ...N spec , 
but also on the particular random distribution of the 
fo.i within the network. Figure 11 displays the differ- 
ence (R) — Rq as a function of the temperature. A log- 
log representation is adopted here for convenience. The 
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dashed straight line represents the best-fit of the data 
with the expression: (R) = R + AT 1 , with 7 = 0.83, 
A = 1.35 x 1CT 2 tl/K~< and R = 0.52 fi. We can see 
that this expression describes rather well the behavior of 
the average resistance. 
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FIG. 11: Variation of the average resistance as a function 
of the temperature. The resistance is espressed in Ohm and 
the temperature in K. The dashed straight line represents the 
best-fit with the expression reported in the figure. 
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FIG. 12: Change of the variance of the resistance fluctuations 
as a function of the temperature. The variance is espressed in 
ft 2 and the temperature in K. The dashed straight line shows 
a best-fit with the expression reported in the figure. 

The change with the temperature of the variance of 
the resistance fluctuations is reported in Fig. 12. Here 
((Ai?) 2 )o is the limit value of the variance at T = 0. 
This value mainly depends on the network size (45| and 
it vanishes for N — > 00. Again, a log- log representation is 
used in this figure. The dashed straight line displays the 
best-fit with the expression: ((Ai?) 2 ) = ({AR) 2 ) + BT r > 
which is found to provide a good description of the data 
for n = 3.5, B = 7.32 x 10" 12 fl 2 /K^ and ((Ai?) 2 ) = 
5.0 x 10" 7 fl 2 . 
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FIG. 13: Relative variance of the resistance fluctuations as 
a function of the average resistance (this last is expressed in 
Ohm). The dashed line shows a best-fit with a power-law of 
exponent s = 2.6. 



Finally, Fig. 13 shows in a log-log plot the relative 
variance of the resistance fluctuations as a function of the 
average resistance. In this case the dashed line represents 
the best-fit with a power-law of exponent 2.6. In other 
terms, Fig. 13 shows that: 



((Ai?) 2 )/<i?) 2 ~ (RY 



(10) 



where s = 2.6. To be consistent with the behaviors of the 
average resistance and of the variance of the resistance 
fluctuations reported in Figs. 11-12, the following scal- 
ing relation should hold within the exponents s, 77 and 7: 
s = 77/7 — 2. By using the previously reported values of 77 
and 7 we obtain for s the value 2.4, in agreement (within 
the error) with the result s = 2.6 obtained by directly 
considering the dependence of ((Ai?) 2 ) / (R) 2 on (R). We 
underline that both, the behavior of the relative variance 
and the value of s, are perfectly consistent with the re- 
sults reported in Ref . [45j . This agreement confirms and 
points out the fact that the model discussed here directly 
generalizes the results of Ref. [45[ to resistors character- 
ized by 1// Q resistance noise. 

IV. CONCLUSIONS 

We have developed a stochastic model to investigate 
the l/f a resistance noise in disordered materials. More 
precisely, we have considered the resistance fluctuations 
of a thin resistor with granular structure in a non- 
equilibrium stationary state. This system has been mod- 
eled as a two-dimensional network made by different 
species of elementary resistors. The steady state of this 
multi-species network is determined by the competition 
among different thermally activated and stochastic pro- 
cesses of breaking and recovery of the elementary resis- 
tors. The network properties have then been studied by 
Monte Carlo simulations as a function of the tempera- 
ture and in the linear regime of the external bias. By 
a suitable choice of the values of the parameters, the 
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model gives rise to resistance fluctuations with different 
power spectra, thus providing a unified approach to the 
study of materials exhibiting either Lorentzian noise or 
l/f a noise. In particular, by increasing the tempera- 
ture in the range from 300 to 800 K the a exponent 
decreases from values near to I down to values around 
0.5, in qualitative agreement with experimental findings 
available from the literature [l|, 0, [nj. This behavior, 
indicative of a drastic reduction of the correlation time 
of the resistance fluctuations, has been also analysed in 
terms of the auto-correlation functions. These last pro- 
gressively change at increasing temperature, going from 



a slow, power-law decay of the correlations at 300 K to a 
faster but non-exponential decay at higher temperatures. 
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